Photonic-crystal slabs with a triangular lattice of triangular holes investigated 

using a guided-mode expansion method 



O 

o 

(N 
S3 

S3 



O 



-a 

S3 
O 

o 



> 

in 

On 

m 
o 

o 
o 



X 



Lucio Claudio Andreani and Dario Gerace* 
Dipartimento di Fisica "A. Volta," Universita di Pavia, via Bassi 6, 1-27100 Pavia, Italy 

(Dated: February 5, 2008) 

According to a recent proposal [S. Takayama et at, Appl. Phys. Lett. 87, 061107 (2005)], the 
triangular lattice of triangular air holes may allow to achieve a complete photonic band gap in 
two-dimensional photonic crystal slabs. In this work we present a systematic theoretical study of 
this photonic lattice in a high-index membrane, and a comparison with the conventional triangular 
lattice of circular holes, by means of the guided-mode expansion method whose detailed formulation 
is described here. Photonic mode dispersion below and above the light line, gap maps, and intrinsic 
diffraction losses of quasi-guided modes are calculated for the periodic lattice as well as for line- and 
point-defects defined therein. The main results are summarized as follows: (i) the triangular lattice 
of triangular holes does indeed have a complete photonic band gap for the fundamental guided 
mode, but the useful region is generally limited by the presence of second-order waveguide modes; 
(ii) the lattice may support the usual photonic band gap for even modes (quasi- TE polarization) 
and several band gaps for odd modes (quasi-TM polarization), which could be tuned in order to 
achieve doubly-resonant frequency conversion between an even mode at the fundamental frequency 
and an odd mode at the second-harmonic frequency; (iii) diffraction losses of quasi-guided modes 
in the triangular lattices with circular and triangular holes, and in line-defect waveguides or point- 
defect cavities based on these geometries, are comparable. The results point to the interest of the 
triangular lattice of triangular holes for nonlinear optics, and show the usefulness of the guided- 
mode expansion method for calculating photonic band dispersion and diffraction losses, especially 
for higher-lying photonic modes. 

PACS numbers: 42.70.Qs, 42.82.Et, 78.20.Bh, 42.65.Ky 



I. INTRODUCTION 



Photonic crystals (PhC) embedded in planar waveg- 
uides, also called PhC slabs, are at the heart of cur- 
rent research on photonic crystals^^^ because of 
the possibility of confining light in all spatial direc- 
tions combined with the advantage of a lithographic 
definition of the pattern. The propagation of light in 
these systems can be controlled by a two-dimensional 
(2D) photonic lattice in the waveguide plane, and by 
total internal reflection in the perpendicular direction. 
Waveguides with strong refractive index contrast (like 
the suspended membrane, or air bridge) support truly 
guided modes lying below the light line dispersion in 
the cladding materials i 7 -* 8 ^ ' 10 ! 11 ! 12 ! 13 ? 14 However, most 
of the photonic modes (or all of them, in the case 
of waveguides with weak refractive index contrast) lie 
above the light line and are only quasi-guided, as 

they are subject to intrinsic losses due to out-of-planc 
diffraction^^d^^^^^.^ 

The physical properties of PhC slabs can be signif- 
icantly different from those of the corresponding 2D 
system, for several reasons. First, the 2D photonic 
modes arc subject to confinement in the vertical waveg- 
uide and the resulting blue shift is strongly polarization- 
dependent. Second, the presence of second- and higher- 
order waveguide modes can produce a complicated pat- 
tern of photonic bands, especially for higher- lying states. 
Third, diffraction losses of quasi-guided modes are an in- 
herent feature of PhC slabs which is absent in the ideal 



2D case. For all of these reasons, some basic and well- 
known properties of 2D photonic crystals cannot be easily 
translated to PhC slabs. For example, it is well known 
that the triangular lattice of circular air holes in 2D sup- 
ports a complete photonic band gap for all propagation 
directions and light polarizations at sufficiently large air 
fractions! 1 ' 26 ^ 7 ' 28 Nevertheless, the same lattice realized 
in a high-index membrane does not possess a complete 
band gap, as the odd modes with respect to a horizontal 
mirror plane (often called quasi-TM modes) are subject 
to a strong and non-uniform blue shift which eliminates 
the gapi^i 2 ^ Indeed, most applications of PhC slabs em- 
ploy the even modes (often called quasi-TE), which do 
possess a band gap for all propagation directions. 

It was recently suggested 3 ^ that the triangular lattice 
of triangular air holes in a high-index membrane gives 
rise to a complete photonic gap for both even and odd 
modes. The physical mechanism is the reduction of sym- 
metry of the basis in the unit cell, as compared to the 
hexagonal symmetry of the 2D lattice, giving rise to a 
splitting between the first and the second odd bands at 
the K point of the Brillouin zone: when realized in a 
high-index membrane, this gap overlaps the usual gap 
between the first and the second even bands, giving rise 
to a complete gap for all directions and polarizations. 
The experimental results reported in Ref. |30| support the 
existence of a complete photonic gap. 

A main purpose of this paper is to perform a system- 
atic study of the triangular lattice of triangular holes, 
as compared to the triangular lattice of circular holes, 
both being realized in a high-index dielectric membrane. 



We calculate the photonic band dispersion, gap maps, 
and intrinsic losses of quasi-guided modes for the 2D lat- 
tice. We also treat line-defect waveguides obtained by 
removing a full row of holes, and point cavities consist- 
ing of three missing holes. In addition to a determina- 
tion of photonic gaps for even and odd modes as a func- 
tion of membrane thickness and air fraction, we com- 
pare diffraction losses for the lattices with conventional 
(circular) and reduced symmetry (triangular) holes: this 
comparison is important in order to assess the possible 
usefulness of the reduced-symmetry lattice. Among the 
results, we find interesting prospects of the triangular 
lattice of triangular holes for nonlinear optics, as it may 
allow to achieve doubly-resonant second-harmonic gener- 
ation (SHG) with an even (quasi-TE) fundamental wave 
and an odd (quasi-TM) harmonic wave when line-defect 
waveguides or photonic cavities are introduced. 

The calculations reported in this paper are performed 
with an approach which we name Guided-Modc Expan- 
sion (GME) method. Maxwell equations are treated by 
expanding the magnetic field into the basis of guided 
modes of an effective homogeneous waveguide, and by 
solving the resulting eigenvalue equation numerically. In- 
trinsic losses of quasi-guided modes are obtained by cal- 
culating the coupling to leaky modes of the effective 
waveguide within perturbation theory (i.e., the photonic 
analog of Fermi Golden Rule for quantum-mechanical 
problems). The GME method, although being an ap- 
proximate one (since the basis of guided modes is not 
a complete basis set), has been applied to a variety of 
photonic lattices and has proven to be useful especially 
for obtaining quasi-guided modes and their diffraction 
lossesi 29 i 31 ' 32 i 33 ' 34 i 35 ' 36 It has also been successfully em- 
ployed for the interpretation of optical experiments on 
PhC slabs ^138 1 39 1 4o 1 4i 1 42 1 43 Another purpose of this paper 
is to provide a detailed description of the GME method, 
together with convergence tests and exemplifying appli- 
cations. 

The rest of this work is organized as follows. In Sec- 
tion II we outline the GME method and discuss a few 
convergence tests. Section III contains the results for the 
2D triangular lattice of triangular holes in a high-index 
membrane, namely photonic dispersion, gap maps and 
intrinsic losses. In Sec. IV we present a few results for 
line-defect waveguides and point cavities in the triangular 
lattice of triangular holes. Section V contains a discus- 
sion of the results and of prospective applications of the 
investigated lattice. Technical details, which are needed 
by the reader in order to implement the GME method, 
are given in the Appendices. 



II. GUIDED-MODE EXPANSION METHOD 

As PhC slabs are intermediate between 2D photonic 
crystals and dielectric slab waveguides, it is reasonable to 
describe photonic modes in these systems starting from 
slab waveguide modes and introducing the effect of a di- 
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Figure 1: Schematic picture of the vertical waveguide struc- 
ture for (a) the photonic crystal slab and (b) the effective 
waveguide. The lower and upper claddings (layers 1 and 3) 
are taken to be semi-infinite. Photonic patterning in the plane 
xy must have the same 2D Bravais lattice for the three layers 
in (a), but it can have different bases in the unit cell. The 
change in grey scales in going from (a) to (b) suggests that the 
effective waveguide is characterized by an average dielectric 
constant in each layer. 



electric modulation in the core and cladding layers. This 
is the central idea of the guided-mode expansion method, 
in which PhC slab modes are expanded in the basis of 
guided modes of an effective homogeneous waveguide and 
coupling to radiative modes is taken into account by per- 
turbation theory. In this Section we describe the formal- 
ism for calculating mode dispersion and intrinsic losses, 
for the general case of an asymmetric PhC slab, and dis- 
cuss convergence of the method. 



A. Formalism for photonic dispersion 

The system we are considering is shown in Fig. [lji. It 
consists of a PhC slab made of three layers 1,2,3, each 
of which is homogeneous in the vertical (z) direction and 
is patterned with a 2D photonic lattice in the xy plane. 
The core layer has a thickness d, while the lower and 
upper claddings are taken to be semi-infinite4^ Pattern- 
ing of each layer is characterized by the same 2D Bravais 
lattice, while the basis in the unit cell can be different - 
a typical situation is that the core layer 2 is patterned, 
while layers 1 and 3 are not (i.e., they can be described by 
the same Bravais lattice as layer 2, but with a vanishing 
air fraction). Writing r = (p,z), the dielectric constant 
e(p, z) is piecewise constant in the z direction and takes 
the form £j(p), j = 1,2,3 in each layer j. We assume a 
magnetic permeability /j, = 1. 

As is well known, Maxwell equations for the electric 
and magnetic fields E, H with harmonic time dependence 
can be transformed into the second-order equation for the 
magnetic field 



V x 



-(r) 



V x H 



4*. 



(i) 



with the condition V • H = 0. By expanding the mag- 
netic field in an orthonormal set of basis states as 



H(r)=^c M H M (r) 



(2) 



3 



the orthonormality condition being expressed by 

y"H;(r)-H,(r)dr = 5 M „ (3) 



(4) 



Eq. ([T]) becomes equivalent to a linear eigenvalue prob 
lem 2 

V 

where the matrix TL^y is given by 

T~t/j.u = J (V x H*(r)) • (V x H„(r)) dr. (5) 

This formulation of the electromagnetic problem bears 
strong analogies to the quantum-mechanical treatment 
of electrons, with the hermitian matrix 7i M „ playing the 
role of a quantum Hamiltonian. Once the magnetic field 
of a photonic mode is known, the electric field can be 
obtained from 



E(r) 



ue(r) 



V x H(r) 



(6) 



and the electric field eigenmodes are orthonormal ac- 
cording to 



/ 



e(r)E*(r)-E„(r) dr = S„ 



(7) 



as is well knowr*^ and easily verified. 

In order to define an appropriate basis set H M (r) for 
the expansion ([2]), we introduce an effective slab waveg- 
uide made of three homogeneous layers with dielectric 
constants t\, e 2 , £3, as illustrated in Fig. [lb. We take Cj 
as the spatial average of tj{p) in each layer: 



4 / £j(p)<lp, 

A J cell 



(S) 



where the integral extends over a unit cell of area A. This 
choice for the average dielectric constants is convenient, 
but by no means unique - numerical tests and other issues 
related to the choice of ej- are discussed in Subsec. lHCl In 
any case, we assume that the average dielectric constants 
fulfill the inequalities 



£2 > Ci, 63 , 



(9) 



in order for the effective slab to support a set of guided 
modes. Their explicit form is well known^ and it will be 
summarized here in order to specify our notation, as well 
as for the paper to be self-contained. Let us denote by 
g = gg the 2D wavevector in the xy plane with modulus 
g and unit vector g, by e g = zxg a unit vector perpendic- 
ular to both g and z, and by uj g the frequency of a guided 
mode which satisfies cg/ri2 < Lu g < eg/ max(ni, 7x3). 
Moreover we define the following quantities: 



Xifl = 
la = 
X3g = 



2\ 1/2 



2 - 

9 -ei-^f 



- w 9 2 

£ 2 — - 9 



1/2 



2\ 1/2 
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(10) 

(11) 
(12) 



which represent the real (imaginary) parts of the 
wavevector in the core (upper/lower cladding), respec- 
tively. The guided mode frequencies are found by solving 
the following implicit equations (with the suffix g being 
understood for simplicity): 

<?(xi + X3) cos(qd) + (xiX3 - q 2 ) sin(qd) = (13) 
for transverse electric (TE) polarization, and 



£2 \ £l 63 / V 6163 



sin(gd) = (14) 



for transverse magnetic (TM) polarization. The guided 
modes at a given wavevector g are labelled by the in- 
dex a = 1, 2, . . . and the eigenfrequencies are denoted by 
Ljga . Explicit forms for the electric and magnetic fields of 
the guided modes are given in Appendix fA) as they are 
needed for the calculation of the matrix elements 

The guided modes of the effective waveguide in Fig. [TJs 
depend on a wavevector g which can take any value in 
the 2D plane. However, photonic modes in the PhC slab 
of Fig. QJ, have the form dictated by Bloch-Floquet theo- 
rem and they depend on a wavevector k which is usually 
restricted to the first Brillouin zone of the 2D lattice: 
indeed, the effect of the dielectric modulation e(r) is to 
fold the guided modes of the effective waveguide to the 
first Brillouin zone and to produce photonic bands and 
band gaps. We therefore write g = k + G, where the 
Bloch vector k lies in the first Brillouin zone and G is a 
reciprocal lattice vector. As the basis states for the ex- 
pansion ^ of the magnetic field, we choose the guided 
modes a (r) of the effective waveguide given in Ap- 

pendix [AJ The guided-mode expansion of the magnetic 
field therefore reads^i 



H k (r) = ^c(k + G,a)H«(r). 



(15) 



The basis set is orthonormal according to Eq. ([3]), but 
not complete, since the radiative modes of the effective 
waveguide are not included in the basis set. This ap- 
proximation will be partially lifted in Sec. IIIB1 where 
the effect of radiative modes in determining diffraction 
losses will be taken into account. 

With this choice for the basis set, the general index /x 
can be written as fi = (k + G, a). The matrix elements 
Ti ^ of Eq. ([5]) can be calculated in a straightforward way 
from the field profiles of the guided modes: the resulting 
analytic expressions are somewhat lengthy and are given 
in Appendix [Bj They depend on the inverse dielectric 
matrices 



Tfe(G,G') 



A 



, \-\ i(G'-G)-p , 



(16) 



in the various layers. Like in usual 2D plane wave calcu- 
lations, a numerically convenient approach to calculate 
these matrix elements is to define the dielectric matrix 



6,(6, G') = - 



i(p)< 



i(G'-G)-p 



dp 



(17) 



and to find rjj(G, G') by numerical matrix inversion as 
T]j(G, G') = e~ 1 (G, G'): this procedure is known to have 



much better convergence properties as a function of the 
number of plane waves ; 48 ' 49 since the truncation rules for 
Fourier series in the presence of discontinuous functions 
are better represented^ If the photonic lattices in layers 
1-3 have a center of inversion (the same for all layers), the 
dielectric matrices £j(G, G ), %(G, G ) are symmetric in 
G, G' and 7i M „ turns out to be a real, symmetric matrix 
when the phases of the fields arc chosen as in Appendix 
lAl In the most general case of a photonic lattice with- 
out a center of inversion (like the triangular lattice with 
triangular holes investigated in this work), the matrix 
Ti-ptu is complex hermitian. Line and point defects can 
be treated by introducing a supercell in one or two direc- 
tions, respectively, like in 2D plane-wave calculations. 

The formalism presented here refers to the most gen- 
eral case of an asymmetric PhC slab. If the structure is 
symmetric under reflection by a horizontal mirror plane 
z = [i.e., if the upper and lower claddings have the 
same dielectric pattern and are made of the same dielec- 
tric materials, ei(p) = e^,(p)\, the basis states as well as 
the photonic eigenmodes can be classified as even or odd 
with respect to the mirror plane. We denote by a xy the 
corresponding specular reflection operator and we shall 
refer to even (odd) states with respect to this operator 
as u xy = +1 (a xy = —1) modes. The simplest way to ex- 
ploit this mirror symmetry is to solve the eigenvalue prob- 
lem Q separately for a xy = +1 and a xy = — 1 states by 
keeping only the appropriate solutions of Eqs. (fl3|) . (fl4|) . 
as detailed in Appendix [Cl There we also discuss the is- 
sue of polarization mixing in the framework of the GME 
method: in general, a photonic mode of a PhC slab is 
a superposition of TE-polarized and TM-polarized basis 
states and all six components of the electric and magnetic 
fields are usually non-vanishing. Another useful symme- 
try property holds when the Bloch vector k lies along 
special symmetry directions, for which mirror reflection 
(Tkz with respect to a vertical plane (fc, z) is a symmetry 
operation of the PhC slab: in this case the eigenmodes 
can be classified as even or odd with respect to this mir- 
ror symmetry, i.e., they have <7kz = +1 or a^ z = — 1 
respectively. Vertical mirror symmetry can be exploited 
by transforming the matrix TL^v with a unitary trans- 
formation which decouples the blocks corresponding to 
(Tkz = ±1 states. 



B. Formalism for intrinsic losses 

When a photonic mode in the PhC slab falls above 
the cladding light line (or light lines, if the waveguide is 
asymmetric), it is coupled to leaky modes of the slab by 
the dielectric modulation and it becomes quasi-guided, 
i.e., it is subject to intrinsic losses due to scattering out of 
the plane. The losses can be represented by an imaginary 
part of the frequency Im(w), which is related to the Q- 
factor of a resonance by Q = w/(2Im(w)). 

Within the GME method, the imaginary part of fre- 
quency can be calculated by time-dependent pcrturba- 
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Figure 2: Schematic picture of radiative (scattering) states 
for the effective waveguide with a single outgoing component. 
States in (a) are outgoing in the lower cladding (layer 1), while 
states in (b) are outgoing in the upper cladding (layer 3), as 
indicated by the thick arrow. The notation for the coefficients 
is appropriate for TE polarization, while for TM polarization 
the replacements W — ► Y, X — > Z have to be made. 



tion theory, like in Fermi Golden Rule for quantum me- 
chanics. The imaginary part of the squared frequency of 
a PhC mode with Bloch vector k, whose frequency lies 
above the cladding light lines (or at least above one of 
them), is given by 



-Im 



= 'E E E 

G' A=TE,TMj=l,3 



/>., 



G' 



(18) 

where the matrix element between a guided and a leaky 
PhC slab mode is given by 



V x Hk(r)) 



(VxHS G , Aj (r))dr (19) 



and pj(k+G'; uif/c 2 ) is the ID photonic density of states 
(DOS) at fixed in-plane wave vector for radiation states 
that are outgoing in medium j: 



P.i 



dk z 

2tt 



-1/2 i 
€■' C 1 



g 2 + kl 



1/2 



4tt 



/ o c 2 q 2 



1/2 



(20) 



Notice the sum over reciprocal lattice vectors and po- 
larizations in Eq. (fT8|) . as all diffraction processes con- 
tribute to Im(w 2 /c 2 ). Equations (fT5|) - ([2U)) generalize the 
expressions given in Ref. [l8| to the case of an asymmetric 
PhC slab and to situations in which processes with all 
G' 7^ contribute to diffraction losses. Once Im(w 2 /c 2 ) 
is found, the imaginary part of frequency is easily ob- 
tained as Im(w) = lm(aj 2 ) / (2oj) . 

For a given wavevector g and polarization A of the ra- 
diative modes, there are two scattering channels, corre- 
sponding to states with an outgoing component in the 
lower cladding (medium 1) or in the upper cladding 
(medium 3)i 45 ' 51 The radiation modes are schematically 
shown in Fig. [2] The state which is outgoing in medium 
1 corresponds to a photonic DOS with j = 1 in Eq. (|20[) . 
while the outgoing state in medium 3 corresponds to a 
DOS with j = 3. Radiation states are normalized ac- 
cording to Eq. ([3]), and only the field components in the 



5 



cladding regions are relevant in determining the normal- 
ization.— Since the field profile of a scattering state tends 
to a plane-wave form in the far field, the photonic den- 
sity of states (|20| is appropriate whatever the explicit 
field profile of radiative PhC slab modes. In the present 
GME method, we evaluate the matrix element by approx- 
imating the radiation modes of the PhC slab with those 
of the effective waveguide. This approximation is consis- 
tent with the treatment of the previous Subsection, as the 
set of guided+radiation modes of the effective waveguide 
is orthonormal according to Eq. ([3]), and will be shown 
later to yield accurate results for diffraction losses when 
compared to more exact methods. In Appendix [D] wc 
give the explicit form of radiation states of the effective 
waveguide which arc outgoing in cither the lower cladding 
(Fig. [2^) or in the upper cladding (Fig.[2jD), for both TE 
and TM polarizations. Using expansion (fTSjl for the mag- 
netic field of a PhC slab mode, the matrix element (fl^|) 
becomes 

"Wk,rad = C(k + G, a)* "Wguided,rad, (21) 

G,a 



where 

H g uidcd, r ad = J ^(Vx[H^ d a (r)*])-(VxHL a + d G , Aj (r))dr 

(22) 

are matrix elements between guided and radiation modes 
of the effective waveguide, which can be calculated ana- 
lytically. In Appendix [E] we collect the matrix elements 
P2"]) . for all possible combinations of polarizations. Those 
expressions are implemented together with Eqs. (fTSj) - 
(|22[) , in order to evaluate the imaginary part of frequency 
and hence the intrinsic diffraction losses of quasi-guided 
modes. 

The calculation of losses is computationally more de- 
manding than the photonic dispersion alone, as the eigen- 
vectors of the linear problem (TJ| are also needed. The 
formalism for intrinsic diffraction losses presented in this 
Subsection and in the Appendix refers to the general case 
of an asymmetric PhC slab. If the slab is symmetric, 
and if the PhC slab modes are found as a xy = ±1 states 
as outlined in Appendix [Cl it is possible to define and 
use symmetric or antisymmetric outgoing states instead 
of the scattering states illustrated in Fig. [2~1 If this is 
done, a xy = +1 (a xy = — 1) quasi-guided modes are cou- 
pled only to symmetric (antisymmetric) radiative states. 
From a computational point of view there is no special 
advantage in using symmetry-adapted radiative states, 
as the computing time needed to implement the pcrtur- 
bative formula Eq. (fTS)) is usually small. Thus even in 
the case of a symmetric PhC slab it is convenient to use 
the general loss formalism, and simply obtain <j xy = +1 
or a xy = — 1 PhC slab states by selecting the basis states 
of the effective waveguide, as explained in Appendix [Cl 



C. Discussion of the method 

In this Subsection we discuss the convergence proper- 
ties of the GME as a function of numerical parameters 
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Figure 3: (a) Real part and (b) imaginary part of frequency as 
a function of number of guided modes N a in the expansion, for 
Cxy = +1 modes in the triangular lattice of circular air holes 
defined in a high-index suspended membrane. The frequencies 
in (a) are calculated at the symmetry points F,K,M, while the 
imaginary parts in (b) are calculated at k — ira/3 along the T- 
K direction. Parameters are: r/a — 0.3, d/a = 0.5, e = 12.11. 
The effective core dielectric constant is fixed to £2 = 8.4827, 
as given by Eq. iJHJ. 



and give a brief overview of previously published results 
and comparison with other approaches. Truncation pa- 
rameters in the GME method are the choice and number 
iVp-w of plane waves and the number N a of guided modes 
of the effective waveguide kept in the expansion. The 
dimension of the eigenvalue problem (TJ} is Apw ■ N a , 
or about half of this value when vertical parity <7k z is 
used. The choice of plane waves can be made in anal- 
ogy to usual 2D plane-wave expansion: the important 
point is how to calculate the inverse dielectric matrix 
r]j(G, G) properly, as already noticed in Subsec. [II Al 
In the present implementation we take an isotropic cutoff 
for reciprocal lattice vectors and adopt the common pro- 
cedure of inverting the dielectric matrix numerically^ 
while this is adequate for many purposes, improved im- 
plementations of the plane- wave expansio n 54 i 55 could also 
be introduced in the GME method in order to speed up 
convergence. 

The optimal N a depends on waveguide parameters and 
on the frequency region considered: typically, for PhC 
slabs defined in a thin high-index membrane, N a = 4 in 
each parity sector (i.e., for both a xy — +1 and a xy = — 1 
modes) yields very accurate results for the lower-lying 
modes. As an example, in Fig. [3] we show the real 
and imaginary parts of the frequency for a few photonic 
modes as a function of the number of guided modes N a , 
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Figure 4: (a) Real part and (b) imaginary part of frequency 
as a function of core dielectric constant in the effective waveg- 
uide, for a xy = +1 modes in the triangular lattice of circu- 
lar air holes defined in a high-index suspended membrane. 
The frequencies in (a) are calculated at the symmetry points 
r,K,M, while the imaginary parts in (b) are calculated at 
k = Tva/3 along the T-K direction. Parameters are: r/a = 0.3, 
d/a — 0.5, e = 12.11. The number of guided modes in the 
expansion is fixed to N a = 4. The vertical bar indicates the 
average dielectric constant, Eq. ([8}, usually employed in the 
GME method (in the present case, £2 = 8.4827). 



for a xy = +1 modes in the triangular lattice of circu- 
lar air holes. Structure parameters are similar to those 
used later in the paper. The number of plane waves is 
truncated to iVpw = 109, which is sufficient for conver- 
gence in the present lattice. It can be seen that both 
the mode energies and the imaginary parts arc stable for 
N a > 3, while a good approximation is obtained already 
for N a = 2 or even N a = 1. This is very convenient 
for the purpose of design and search of parameters, as a 
photonic structure in a high-index contrast PhC slab can 
be optimized within a fully 3D approach, with the same 
computing time of a conventional 2D plane-wave expan- 
sion. The precise calculation (e.g., with N a > 4) can be 
performed a posteriori once the structure parameters are 
optimized. 

It is also important to discuss the choice of the dielec- 
tric constant of the effective waveguide, which is usually 
taken to be the spatial average |(5J) in each layer. While 
this choice is a natural one, it might be asked whether 
the results depend on the assumed value of Ej. In Fig.Q] 
we show the real and imaginary parts of frequency for 
a few modes, with the same parameters of the previous 
Figure, as a function of the effective dielectric constant in 
the core region. It can be seen that the results are rather 
insensitive to the assumed value of (=2, except when it be- 



comes smaller than about five. All curves are flat around 
the value li = 8.4827 which corresponds to the spatial 
average. This satisfactory behavior arises because the 
GME method employs the set of all guided modes of 
the effective waveguide (although in practice they are 
truncated to a finite number for numerical convenience) , 
thereby compensating for different choices of the effective 
waveguide used to define the basis set. Thus the choice 
of the dielectric constant of the effective waveguide as 
the spatial average of e,(r) is justified, at least for the 
low air fractions that are commonly employed. Although 
one might consider extending the method by introduc- 
ing an anisotropic dielectric tensor for each layer of the 
effective waveguide, Fig. @] strongly suggests that the re- 
sults of the GME method would not be changed by such 
a complicated extension which is therefore unnecessary. 

The approximations of the GME method are: (i) for 
photonic dispersion, the shift of guided and quasi-guided 
modes due to second-order coupling to leaky modes of the 
effective waveguide is neglected; (ii) for diffraction losses, 
the density of radiative modes of the PhC slab is approx- 
imated with that of the effective waveguide. The effect 
of these approximations can be judged by comparing the 
results with those of other methods that are known to be 
exact within numerical accuracy. A full analysis of this 
issue is outside the scope of this paper, however com- 
parisons between the results of GME and those obtained 
by other methods have already been published and are 
shortly summarized here. The dispersion of quasi-guided 
modes was found to be in very good agreement with the 
resonance positions in reflectance or transmittance spec- 
tra found from a scattering-matrix treatment^ In fact, 
the complex frequencies of photonic modes can be ob- 
tained from the poles of a scattering matrix i 19 i 22 i 23 The 
intrinsic propagation losses of line-defect waveguides in 
the triangular lattice obtained with different methods 
(GME^ finite-difference time-domain or FDTD 2 ^ and 
from the poles of a scattering matri x 19 ' 56 ) were compared 
for the same waveguide parameters and were found to 
agree very well with each other: a detailed comparison is 
presented in Ref. [H. A comparative study of GME and 
Fourier modal methods applied to cavity modes in one- 
dimensional PhC slabs has been reported in Ref. H^, and 
very good agreement was found for the photonic disper- 
sion as well as the intrinsic losses. Finally, the Q-factors 
of cavity modes in LI, L2 and L3 nanocavities (one, two 
or three missing holes in the triangular lattice) as a func- 
tion of nearby hole displacement arc found to be almost 
identical when calculated with either the GME method^ 5 - 
or the poles of a scattering matrix ^ Agreement with 
experiments performed on PhC slabs made of Silicon- 
on-Insulator ,22^^ Si membranes^ GaAs/AlGaAs^^ 
and Silicon Nitride^ is also very satisfactory. The set of 
these comparisons indicates that the GME is a reliable 
method for both photonic mode dispersion and intrin- 
sic losses, whose main advantages are computational effi- 
ciency and ease of application to various kinds of vertical 
waveguides and of photonic lattices, both periodic in 2D 
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Figure 5: Schematic picture of the triangular lattice of trian- 
gular air holes with lattice constant a and triangle side L; the 
unit cell in real space is indicated. The lattice is patterned in 
a high index suspended membrane of thickness d. The Bril- 
louin zone with main symmetry directions in reciprocal space 
is also drawn. 

and containing line- and point defects. 

III. TRIANGULAR LATTICE WITH 
TRIANGULAR HOLES 

We report in this Section a systematic study of the tri- 
angular lattice of triangular air holes patterned on a free- 
standing high-index slab of thickness d, as schematically 
shown in Fig. [5l The relevant dimensions of this lattice 
are the triangular basis side, L/a, and the slab thick- 
ness. Throughout this Section we assume (unless other- 
wise specified) that the lattice is patterned in a material 
with dielectric constant e = 12.11, which is appropriate 
for Silicon at the usual telecom wavelength A = 1.55 /im 
but also to other materials with similar dielectric con- 
stants (e.g., GaAs). For the calculations shown in the 
present Section, we employed up to 109 plane waves and 
4 guided modes per parity sector in the basis set of the 
GME method. All the calculations have been done on 
an ordinary Pentium computer with low computational 
effort. The complex Fourier transform of the dielectric 
constant for the triangular lattice of triangular holes is 
given in Appendix [F] Notice that the lattice is invari- 
ant under rotations by 120° and is not centrosymmctric, 
but since u>(— k) = w(k) by time-reversal symmetry, the 
same symmetry directions of the usual triangular lattice 
of circular air holes do apply. 

The first realization of such a lattice has been reported 
in Ref. l3(il. with the purpose of measuring a complete pho- 
tonic band gap in a waveguide-embedded 2D PhC thanks 
to a reduction of symmetry for the basis in the unit 
cell of the triangular lattice, as already pointed out in 
the Introduction. The importance of finding a complete 
photonic band gap in a waveguide-embedded photonic 
structure is related to the practical realization of linear 
waveguides or cavities with reduced losses or higher Q- 
factors, as scattering into quasi-guided modes of opposite 
parity would be avoided. In Fig. [6] we present the pho- 
tonic band dispersion calculated with the GME method 
for structure parameters as in Ref. [3(3. Even (TE-like) 
and odd (TM-likc) modes are reported in Figs. [5^ and 




Figure 6: Photonic band dispersion for (a) even {a xy = +1) 
and (b) odd {a xy = —1) modes of the PhC slab with a tri- 
angular lattice of triangular holes, parameters as in Ref. |3C| : 
L/a = 0.85, d/a = 0.68, e = 12.2. Light dispersion in air is 
represented by dotted lines. 



b, respectively. It is worth noticing that in the absence 
of second-order waveguide modes, the fundamental gaps 
of both parities would overlap in the dimcnsionlcss fre- 
quency range 0.31-0.32, as it can be seen from the Fig- 
ure. However, the presence of a second-order mode with 
cut-off around uia/(2iTc) ~ 0.28 for both polarizations, 
as shown in Fig. [5J prevents the opening of a true pho- 
tonic band gap. Such second-order modes (which do not 
appear in a 2D treatment with an effective index) are un- 
avoidable in the present structure, due to the relatively 
large slab thickness, d/a ~ 0.7. 

We believe that our result is not in contrast with the 
experimental findings reported in Ref. [30l . in which a 
strong reduction of the transmission intensity was mea- 
sured along the TK and TM directions. Since the disper- 
sion of the second-order waveguide mode is very close to 
the air light line in the relevant frequency region, such 
modes are barely coupled to the external beam focused 
on the sample side and also they are subject to large 
disorder-induced losses; thus, it is likely that transmis- 
sion into such modes is attenuated in a few lattice pe- 
riods, which would explain the apparent photonic band 
gap observed experimentally. On the other hand, when 
line- or point-defects are introduced, such second-order 
modes are necessarily present within the gap, thus lim- 
iting the usefulness of this design for applications. In 
the following, we will consider a photonic band gap for 
the waveguide-embedded structure only when no higher- 
order modes are present within the apparent gaps formed 
by the fundamental modes, i.e., when the PhC slab is 
truly single-mode. 

With the aim of optimizing the design proposed by 
Takayama et at, we performed a systematic study of the 
present photonic lattice at varying structure parameters. 
In Fig. [7] we show the photonic band dispersion and gap 
maps for the triangular lattice of triangular air holes with 
slab thickness d/a = 0.3 (Figs. [7^ and d), d/a = 0.5 
(Figs. [TJj and e) , and for the reference 2D system of in- 
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L/a Ua Ua 

Figure 7: (Color online) Upper panels: photonic band dispersion for L/a = 0.8 in (a) PhC slab of thickness d/a = 0.3, (b) 
d/a — 0.5, and (c) ideal 2D photonic crystal. Lower panels: gap maps as a function of the holes' side for (d) d/a = 0.3, (e) 
d/a = 0.5, and (f) ideal 2D system. The high index material is assumed in all of these calculations to have a dielectric constant 
e = 12.11. Light grey regions represent the modes above the light line for the waveguide-embedded systems. 



finite thickness (Figs.[7j; and f), respectively. As is well 
known for more conventional lattices^^ the blue shift 
of photonic modes with respect to the ideal 2D system 
due to the dielectric confinement is strongly polarization- 
dependent; in particular, it is more pronounced for odd 
(TM-like) modes. This physical behavior is seen by com- 
paring Figs. [JJi-c, in which the dispersion of photonic 
modes is calculated for the same 2D photonic lattice in 
which triangular air holes have side L/a = 0.8. For 
d/a = 0.3 the second-order cut-off is rather high in en- 
ergy, thus the system can be considered as single-mode 
below the light line. On increasing the slab thickness, 
the second-order cut-off red shifts. For d/a = 0.5, both 
even and odd modes open their fundamental band gaps, 
but they do not overlap to give a complete photonic 
band gap for all symmetry directions. In order to span 
a wider range of structure parameters, important infor- 
mation can be inferred by the gap width as a function 
of the holes' side. As it can be seen in Figs. [7Ji-f, the 
gap maps depend strongly on the waveguide thickness 
and differ substantially from those of the 2D system. For 
d/a = 0.3, the fundamental odd gap is closed for all val- 
ues of L/a, due to the larger blue shift of band edges at 
the K-point than at the M-point. Even if both funda- 
mental band gaps are present for d/a = 0.5, no complete 
gap opens due to the stronger blue shift of odd modes. 
For d/a > 0.5, the presence of higher-order modes pre- 
vents the opening of complete gaps as well. The gap 
map for the ideal 2D system presents a great variety of 
photonic gaps for both polarizations, and also some com- 
plete gaps at high frequencies. It is worth noting that 
the fundamental TE and TM band gaps do not overlap 
in frequency, unlike in the ordinary triangular lattice of 



circular air holesi Furthermore, the fundamental TE gap 
turns out to be narrower than the usual gap for circular 
holes. 

Even if no truly complete band gaps appear to be 
present for this waveguide-embedded lattice, it is in- 
teresting to investigate the formation of resonant gaps 
at lu and 2u>, to be exploited in non- linear optical 
applications like second- harmonic generation (SHG). 
While doubly-resonant SHG has been widely stud- 
ied in one-dimensional systems like multilayers and 
microcavities,^a^^i^^^^ no 2D photonic lattice 
giving such an interesting result has been reported in the 
literature. In particular, it can be seen from Fig. [TJd that 
many odd gaps are present for a slab thickness d/a = 0.5. 
In this case, it is easy to find a condition for which the 
fundamental o~ xy = +1 gap at pump frequency to is res- 
onant with a a xy = — 1 gap at 2u>. One of the possi- 
bilities is indeed shown in Fig. [8] It is substantially the 
same result reported in Fig. [7b, but here the odd modes 
(Fig. [5}j) are plotted in a doubled frequency range with 
respect to the even modes (Fig. [5^)- As it can be seen, 
there is a good overlap in frequency between the large 
band gap of Fig. [5^ and the rather narrow but robust 
gap of Fig. [5Jd. This could be interesting, e.g., for appli- 
cations to non-linear optical converters based on GaAs, in 
which the pump s-polarizcd field would be converted into 
a p-polarized second-harmonic fields In this respect, the 
usefulness of the present lattice would be enhanced with 
ad hoc fabricated line- or point-defects. It should be re- 
minded that these results have been obtained by using 
the same dielectric constant for both pump and second- 
harmonic modes, but it is easy to extend the present cal- 
culations to non-linear materials with optical constants 
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of known frequency dispersion. 

One of the issues when dealing with the triangular lat- 
tice of triangular holes would be given by intrinsic losses. 
In order to show the general interest of the triangular 
hole-based lattice, and also to show an application of the 
GME method, we compare intrinsic losses of the present 
system to the ones of an ordinary circular hole-based PhC 
slab. In this case, we restrict our discussion to <r xy = +1 
modes along the TK symmetry direction. Intrinsic losses 
are quantified by calculating the imaginary part of fre- 
quency, as detailed in the previous Section. In Fig.[9^i and 
c we show the photonic dispersion and intrinsic losses for 
the first few bands of an ordinary PhC slab with a tri- 
angular lattice of circular holes. It can be noticed that 
the modes along the chosen symmetry direction are sepa- 
rated according to their mirror symmetry with respect to 
the vertical plane, which can be cither even {o^z = +1) or 
odd ((Tfer = — On the right hand side, in Fig.O) and 
d, we plot dispersion and losses for the case of triangu- 
lar holes. In this case no symmetry operation <7k 2 holds, 
so the modes are mixed. The choice of r/a = 0.3 for 
the circular holes and L/a = 0.8 for triangular ones, re- 
spectively, allows to obtain roughly the same air fraction 
/ ~ 0.32 for the two lattices. This can be also verified by 
the roughly equal second-order cut-offs. The dispersion is 
only slightly modified when changing the hole shape from 
circular to triangular, and in particular it appears that 
a general blue shift occurs for the modes in Fig. [SJj with 
respect to the ones in [5^. For what concerns intrinsic 
losses, we can see that they remain almost of the same or- 
der on going from circular to triangular holes, at least for 
the lower frequency bands. This means that (in the ab- 
sence of disorder) no particular differences are introduced 
by changing the basis of the lattice, apart from the loss of 
symmetry tJkz ■ As the band labelled with 1 is almost un- 
changed on going from the circular to the triangular hole 
shape, also the intrinsic losses are quite similar between 
the two lattices. For higher lying bands, in particular 
for the band labelled with 2, such similarity becomes less 
evident on approaching the light line, owing to the blue 
shift of modes for the triangular hole-based lattice and 
to the mixing with other modes which does not occur in 
the conventional lattice. In particular, owing to the lack 
of Okz symmetry, there is an anticrossing between mode 
2 and the second-order waveguide mode occurring very 
close to the light line, which correspondingly affects the 
imaginary part of frequency (see Figs. 0>d). The mode 
labelled with 3 has larger losses for triangular holes than 
for circular ones, but it is quite high in energy. 

As a conclusion to this section, we have shown that the 
triangular lattice of triangular air holes does indeed have 
a complete photonic band gap between the fundamental 
waveguide modes, but the gap is actually eliminated by 
the presence of second-order modes. The lattice can be 
designed for nonlinear optical applications in order to 
have a gap for even modes at the pump frequency which 
is resonant with a gap for odd modes at the SH frequency. 
Intrinsic losses arc of the same order as in the circular 




Figure 8: Photonic band dispersion for (a) a xy = +1 modes 
at pump and (b) a xy = — 1 modes at second- harmonic fre- 
quencies, respectively, showing resonant band gaps at uj and 
2uj useful for non-linear applications. Parameters are as in 
Fig. [7b: L/a = 0.8, d/a = 0.5, e = 12.11. Light dispersion in 
air is represented by dotted lines. 

hole-based triangular lattice. The GME method is useful 
for a comparison of the two lattices and for the design of 
structures with desired characteristics. 



IV. LINE-DEFECT WAVEGUIDES AND POINT 
CAVITIES 

We show in the present Section that line- and point- 
defects can be introduced into the triangular-hole lat- 
tice without any loss of performance with respect to the 
circular-hole one. It is well known that removing a row of 
holes in a 2D photonic lattice introduces a linear defect 
supporting guided modes within the photonic band gap. 
In a triangular lattice such waveguide is usually called 
Wl, and its channel width (namely the distance between 
the holes surrounding the waveguide) is w = wq = y/3a. 
These guided modes exploit the usual index confinement 
along the vertical direction added to the gap confinement 
in the waveguide plane, and they are studied with great 
interest by the scientific community owing to prospective 
applications in photonic integrated circuits. 

In Fig. [TU] wc report a comparison between the Wl 
waveguide realized in the triangular lattice with circular 
fFigs. ITUh-c) and triangular holes (Figs. [TUb-dh respec- 
tively. For the circular hole-lattice the guided modes are 
classified according to their mirror symmetry with re- 
spect to the vertical plane kz. For such calculations, we 
used supcrcells along the TM direction of widths rang- 
ing from 4u>o to 9wq, with 181 to 323 plane waves and 4 
guided modes in the basis set for the expansion. Then, 
an average was taken to get final results for losses, in or- 
der to smooth out finite supercell effects^ As it is seen 
from the Figure, the dispersion of the guided modes is 
very similar between the two systems, apart in the re- 
gion of mixing. In fact, as o^z is a symmetry operation 
for the system in Fig.fTUh. the two guided modes are com- 
pletely decoupled, while in Fig. 110b they mix producing 
a consistent anti-crossing and a mini-gap. Looking at 
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Figure 9: (a) Photonic band dispersion and (c) intrinsic losses 
for tj xy = +1 modes in the triangular lattice of circular air 
holes with r/a — 0.3, d/a — 0.5, e = 12.11, along the T-K 
symmetry direction; modes are classified according to their 
Bkz parity, (b) Dispersion and (d) intrinsic losses for the 
same lattice of triangular air holes with L/a = 0.8, d/a = 0.5, 
e = 12.11. Light dispersion in air is given by dotted lines. 



the imaginary part of frequency, we see by comparing 
Figs. [TUb and d that they are almost identical away from 
the mixing region, which means that Wl waveguides in 
triangular lattice of triangular holes may behave as stan- 
dard Wl waveguides j££ Such result is quite surprising, as 
the asymmetry of such a waveguide could be expected to 
lead to higher scattering of propagating light states. On 
the contrary, the reduction of symmetry does not increase 
losses in the perfect photonic lattice - of course, the role 
of fabrication imperfections may affect differently the two 
structures in real lattices. 

The final goal of this paper is to exploit the potential of 
the triangular lattice of triangular holes for non-linear op- 
tical applications. In this respect, the ultimate system for 
frequency conversion at the nanoscale is a photonic cav- 
ity with double resonance for pump and harmonic waves. 
Doubly-resonant microcavities for second-harmonic gen- 
eration have been studied in the context of Fabry-Perot 
cavitiesi 58 i 59 ' 63 i 64 ' 65 According to a recent proposal, dou- 
ble resonance in photonic crystal cavities may allow to 
achieve strong coupling between single photons^ It has 
been shown that full 3D confinement of light with a very 
high figure of merit (quality factor, Q) can be achieved 
by properly designed point-defects in 2D PhC slabs' 
In particular, removing three holes along the TK direc- 
tion (L3 cavity) and optimizing the positions of the holes 
surrounding the cavity allows to achieve Q-factors of the 
order of lO 5 ^ We compare here the L3 cavity structure 
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Figure 10: (a) Photonic band dispersion and (c) intrinsic 
losses for a xy = +1 modes in Wl waveguides in triangu- 
lar lattice of circular air holes with r/a — 0.3, d/a — 0.5, 
e = 12.11; modes are classified according to their atz parity, 
(b) Dispersion and (d) intrinsic losses for the same waveguide 
in the triangular lattice of triangular air holes with L/a = 0.8, 
d/a = 0.5, e = 12.11. Light dispersion in air is given by dot- 
ted lines. The shaded regions represent the bulk modes of 
the respective lattices projected onto the line-defect Brillouin 
zone. 



made in a triangular lattice of triangular air holes to the 
usual L3 cavity with circular holes and the results are 
shown in Fig. 1111 For such calculations, a square su- 
percell of dimensions 10a x 5wo has been taken. Up to 
1551 plane waves and two guided modes were used in 
the GME. As it can be seen from Fig. [TTa . many dis- 
pcrsionless <j xy = +1 defect modes appear in the pho- 
tonic band gap, due to the tight confinement provided 
by the cavity. We can evaluate the vertical Q-factor of 
such modes (i.e., the one determined by intrinsic out-of- 
plane losses) by the definition Q = oj/[2lm(ui)], where 
both lu and Im(u>) have been averaged over the folded 
Brillouin zone in order to reduce the effects of a finite 
supercell size. Very good agreement with measured Q- 
factors and defect mode energies has been already found 
with this GME-based procedure^ The Q-factor of the 
un-optimized L3 cavity in the lattice of circular holes is 
theoretically estimated to be Q ~ 5400. When the same 
cavity is realized in the triangular lattice of triangular 
holes (Fig. [Tib), we see that the fundamental mode is 
almost unchanged, and its Q-factor is comparable to the 
one of Fig. [TTa . The slight blue shift of the defect mode 
and of the lower band edge in Fig. [TTb with respect to 
Fig. [TTk is due to the slightly different air fractions of 
the two lattices. As we have already shown in Fig. [5J 
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Figure 11: Cavity modes for (a) L3 cavity in triangular lattice of circular holes, a xy = +1 modes, with r/a — 0.3, d/a — 0.5, 
e = 12.11; (b) L3 cavity in triangular lattice of triangular holes, a xy = +1 modes at pump frequency, with L/a = 0.8, d/a — 0.5, 
e = 12.11; (c) L3 cavity in triangular lattice of triangular holes, a xy = —1 modes at second harmonic frequency. Shaded regions 
represent the bulk modes of the respective lattices projected onto the point-defect Brillouin zone. 



we can design a triangular lattice with triangular holes 
having a doubly-resonant band gap. Thus, we plot in 
Fig. [TTb the dispersion of a xy = — 1 modes within the 
doubled frequency range with respect to Fig. lllb . in anal- 
ogy with Fig.[8j Many defect modes appear within the 2uj 
band gap, and the lowest mode has a Q-factor of almost 
10 3 , which is a very promising value considering that the 
structure has not been optimized in any way. There is, 
in fact, much room for improving the Q-factors at both cu 
and 2lo by means of geometry optimization of the nearby 
holes. The present result make us confident that high- 
Q, doubly-resonant nanocavities may be realized in PhC 
slabs made of non-linear materials, bringing new exciting 
results in integrated optics at nanoscale level. 



V. CONCLUSIONS 

The triangular lattice of triangular holes in a high- 
index dielectric slab has been thoroughly investigated 
by means of the guided-mode expansion method, whose 
detailed formulation and examples of applications have 
been given in the paper. The investigated lattice does 
have a complete band gap for both polarizations when 
the fundamental waveguide mode is considered, consis- 
tently with the results of Ref. [3(| but the gap is actu- 
ally eliminated by the presence of second-order waveguide 
modes. The intrinsic losses of quasi-guided modes in the 
triangular lattice with triangular and circular holes, and 
in line-defect waveguides as well as nanocavities defined 
therein, arc comparable. 

An interesting application of the triangular lattice of 
triangular holes for nonlinear optics follows from the fact 
that a gap for a xy — +1 (TE-like) modes at the pump 
frequency can be designed to be in resonance with a 
gap for a xy = — 1 (TM-like) modes at the harmonic fre- 
quency: thus, resonant second-harmonic generation ex- 



ploiting band-edge effects may be realized. Promising re- 
sults are also found for point cavities: the L3 cavity in the 
triangular lattice of triangular holes supports a u xy = +1 
cavity mode at to which can be resonant with a a xy = — 1 
cavity mode at 2lo. There is much room for improving 
the design of doubly- resonant, high-Q nanocavities by 
geometry optimization. These results are interesting in 
view of application of the triangular lattice of triangular 
holes for nanoscale nonlinear optical processes. 
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Appendix A: GUIDED MODES OF THE 
EFFECTIVE WAVEGUIDE 

The frequencies of the guided modes of the effective 
waveguide are found from the implicit equations Q13p 
for TE polarization, and (fH)) for TM polarization. The 
guided modes are labelled by the wavevector g and the 
mode index which can be combined into a single 
index /j, = (g, a). The mode profiles for TE polariza- 
tion can be written as follows (with the time dependence 
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being understood), with reference to the geometry of Fig. [T] 
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In the above expressions, S is a normalization surface, 
which disappears from the final results, while the mag- 
netic field is obtained from the electric field thorugh the 
Maxwell equation H(r) = —77V x E(r) . Application 
of standard transfer-matrix theory leads to the follow- 
ing relations between coefficients (with the suffix \i being 
understood for simplicity): 
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and to the well known implicit equation (|13|) of the main 
text. The normalization integral Q implies the relation 
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which, together with Eqs. (|A3[) - (|A5[) , determines all co- 
efficients Aj and Bj. 

The guided modes for TM polarization are character- 
ized by the following field profiles: 
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(p> 2 ) = -7=-*— 



C 2 ^ z + D 2 „e~ 
jDlMe xi M ( 2 +d/2) i 



"hi - 



Z > 
\z\ < 

z < - 



^Cs^ X 3„g + igB)e-^- d/2 \ 

X[CV(-g M § + gz)e lq " z + D 2 „i{q„g + gz) e - lq * z ], 
^D^(- Xlll g + igz)e^ +d / 2 \ 

I 



z > 
\z\ < 



z < 



(A7) 
(A8) 



where the electric field follows from the magnetic field 
profile through Eq. The appropriate relations be- 

tween coefficients can be derived from (|A3I) - (|A5|) by the 
replacements A — ► C, B — ► D and xi - * Xi/^ii Q ~ * l/^2 
(except in the trigonometric functions), X3 — * X3/^3- The 
normalization integral of the magnetic field gives the con- 
dition 



\Di\ 



2xi 2 X3 



+ J 



-+d 



\c 2 \ 2 + \d 2 \ 2 + (c;d 2 + d* 2 c 2 ) 



sin(gd) 



qd 



(A9) 



Appendix B: MATRIX ELEMENTS FOR 
PHOTONIC DISPERSION 

We calculate the matrix elements Ti^ of Eq. ([5]) 
between guided modes of the effective waveguide from 
the field profiles given in Appendix |A] for TE and TM 
polarizations. We adopt the following notations: /1 = 



(k + G, a) = (g, a) and v = (k + G', a') = (g', a'), where 
k is the Bloch-Floquet vector, G and G' are reciprocal 
lattice vectors and a, a' are the indices of guided modes 
of the effective waveguide at wavevectors k + G, k + G' 
respectively. There are four kinds of matrix elements: 
TE-TE, TM-TM, TE-TM and TM-TE. The z integral in 
Eq. ([5]) can be broken into three terms over the regions 
1,2,3 where the dielectric constant does not depend on z, 
and can be expressed in terms of the following integrals: 



h 
h 



d I 2 
d/2 



poo 

/ e -(«.+X3,)(-d/2) d2 = fefi + Xa J-l t 
Jd/ 2 

L 
I 



i(( qfl ±q v )d/2) 



d/2 
-d/2 



3 ( ^ +x - )(z+d/2) d 2 = (xiM+xi,r 1 . 



The xy integrals in Eq. ([5]) yield the inverse dielectric 
matrices in the three layers, T)j(G, G ), defined in Eq. 
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[[TBI). The matrix elements arc then obtained as follows: 



Hi 



+{l 3 ) 2 m (G,G')A^A z „h + {l 2 fr, 2 {G,G') 
x [(A 2fi A 2u + B2„B 2 „)I 2 - + (A 2 „B 2v , + B 2 ^A 2v )I 2 +\ } , 

(Bl) 



H 



= rji (G, G')Dl il Dx„(xiij.Xiv9 • 9 + 99 )h 



+m{G, G')C3 ll C 3v (x3nX3i>g ■ g + gg')h 

+ V2 (G,G') [(C^C 2v +D* 2 ^D 2 „) 

x {Wv9 ■ 9 + gg')h- + (C 2fJ ,D 2v + D* 2 ^C 2v ) 

x {-g^g ■ g + gg')h+\ , (B2) 



+e :i rj :i (G, G')A* Zll C^XAvh + ie 2 r] 2 (G, G')q„ 
x [(-A^C 2v + B^D 2v )h- 
+(A^D 2u - B* 2 ^C 2v )I 2+ }} , 



(B3) 



= y—) 9 -e g ' {-ei»?i(G, G^-D^-B^xiaJi 

+£37?3(G, G')C3 m A3„X3m^3 

-il 2 r) 2 (G, G')q„ [{-C 2li A 2v + D 2fi B 2l/ )I 2 - 

+ {D* 2 ^A 2v - C^B 2u )I 2+ ) } . (B4) 



Appendix C: SYMMETRY PROPERTIES AND 
TE/TM MIXING 

If the PhC slab is invariant under reflection through a 
mirror plane bisecting the slab (which we take as z = 0) , 
the effective waveguide has the same reflection symmetry 
and its eigenmodes can be classified as even or odd with 
respect to the reflection operator a xy defined in the main 
text. The spatial symmetries of the nonvanishing field 
components are given in Table UJ assuming a wavevector 
along the £ direction. The guided mode frequencies are 
found from the following implicit equations: 



qd 
T 

qd 
i — 
2 



qd 

T 

qd 
~2 



e 2 
9_ 



2 ei 2 

qd xi Qd 

cos 

2 ei 2 



= 


TE, a xy 


= +1 


(CI) 


= 


TE, a xy 


= -1 


(C2) 


= 


TM, a xy 


= +1 


(C3) 


= 


TM, (J X y 


= -1 


(C4) 



Notice that the first-order TE and TM modes of the ef- 
fective waveguide are described by Eqs. (|C1|) and (|C4[) . 
respectively. By keeping only a xy = +1 (a xy = —1) ba- 
sis states in the guided-mode expansion, photonic eigen- 
modes which are even (odd) with respect to horizontal 



Table I: Spatial symmetry of the electric and magnetic field 
components of an eigenmode of a symmetric waveguide with 
respect to the mirror symmetry operation a xy . The wavevec- 
tor is assumed to lie along the x direction. Vanishing compo- 
nents are denoted by "*". 





E x 


Ey 


E z 




Hy 




TE, 


@xy 


= +1 


* 


+ 


* 




* 


+ 


TE, 


@xy 


= -l 


* 




* 


+ 


* 




TM, 


@ xy 


= +l 


+ 


* 




* 




* 


TM, 


&xy 


= -l 




* 


+ 


* 


+ 


* 



mirror symmetry arc obtained, and the computational 
effort in solving the eigenvalue equation ^ is reduced. 

In general, a photonic eigenmode in a PhC slab is 
a linear combination of TE- and TM-polarizcd basis 
states and all six components of the electric and mag- 
netic field arc nonvanishing. Nevertheless, low- lying pho- 
tonic modes are often dominated by the lowest-order 
guided mode of the effective waveguide. In this situation, 
o xy = +1 states are dominated by the TE guided modes 
described by Eq. I|C1[) and can be called "quasi-TE" , 
while <j xy = — 1 states arc dominated by the TM guided 



modes described by Eq. (|C4|) and can be called "quasi- 
TM" . This widespread terminology is useful and appro- 
priate for large wavevectors, when the dominant field 
components are the spatially even ones: making reference 
to TableUJ and since the dielectric modulation €j (xy) cou- 
ples guided modes with all wavevector directions in the 
2D plane, the dominant components are E Xl E y , H z (for 
<j xy — +1 or TE-like modes) or H Xl H y , E z (for a xy = — 1 
or TM-like modes). However, the terminology becomes 
inadequate when the PhC slab is multimode, as it always 
happens at sufficiently high frequencies. If the waveguide 
is not symmetric, any PhC slab mode is a linear combina- 
tion of basis states arising from all equations (|Cl|) - (|C4p 
and even the lowest-order TE/TM guided modes of the 
effective waveguide are mixed. 



Appendix D: RADIATION MODES OF THE 
EFFECTIVE WAVEGUIDE 

In this Appendix we give the radiation modes of the 
effective waveguide which are outgoing in either the lower 
or the upper cladding (see Fig. (2). Each of these radiation 
modes is labelled by the wavevector g in the 2D plane, 
the frequency to, and the polarization. In the following 
the quantum numbers g, to will be understood in order to 
simplify the notations. Instead of distinguishing between 
real and imaginary components of the wavevectors in the 
three regions as in Eq. (|12p . we define three complex 
wavevectors as 



1/2 



1,2,3. 



(Dl) 
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While the wavevector g2 in the core region is always 
real, the wavevectors q\ , q^ in the cladding regions can be 
either real or purely imaginary according to the values 
of the dielectric constants and of g,u>. If the effective 
waveguide is symmetric (ei = €3) all radiation modes 
have real (71,93, while if t\ > €3 there is a frequency 
region cg/^/ei < to < cg/^/ei in which q\ is real but q% 
is imaginary: i.e., a quasi-guided mode in this region of 
k-u) space is diffracted to the lower (but not to the upper) 
cladding. This phenomenon is automatically taken into 
account by the Kronecker ^-function in the DOS formula 
(|20p . as the photonic DOS vanishes when qj is imaginary. 
Therefore we shall write down the electric and magnetic 
field components of radiative waveguide modes assuming 
that all gj's are real - if one of them is purely imaginary, 
that radiation mode does not carry an energy flux and it 
docs not contribute to scattering loss. 

With this proviso, the radiation modes of the effective 
waveguide for TE polarization in the three regions j = 
1, 2, 3 can be written as follows (setting z\ = —d/2, z^ = 
0, Z3 = d/2 and with the time dependence e~ %UJgt being 
understood): 



ETW) 



Igp 



k [ 



iqj (z-Zj) 



], (D2) 



H5"W) = \{gz- qj g)W^ iqj{z - Zj) 
VS u L 



(D3) 



For the scattering state which is outgoing in the lower 
cladding (Fig. [2^,), W\ = and the normalization con- 
dition ([3]) determines the coefficient of the outgoing 
component^ as X\ = l/y/e[. All other coefficients are 
then found from a standard transfer-matrix calculation, 
which yields 



(92 - gi)e- S92d/2 (q 2 + qi)e-^ 2 



' w 2 " 


1 




~ 2q~2 


and 




" w 3 " 


1 









' Wi " 







(D4) 



(q 3 +q 2 )e i «* d / 2 (53 - q a )e- i * d ' a 
{q,-q2)e^ 2 (qs + q2)e-"* d / 2 



W 2 
X 2 
(D5) 

On the other hand, for the state outgoing in the upper 
cladding (Fig. [2Jd), A3 = and the proper normalization 
condition^ is W3 = \/y/€$. The other coefficients can 
then be found from the inverse relations of (|D4|) . (|D5|1 . 

For what concerns TM-polarized radiation modes, the 
field profiles are given by 



H 5 ad (/=>>«) = 



>], (D6) 



VS ejU L 

+ {gz + qfg)Z 3 e-^-^]. (D7) 



For the scattering state outgoing in the lower cladding 
(Fig. [5^), Y\ = and the normalization condition ([3]) 
determines the coefficient of the outgoing component as 
Z\ = 1. For the state outgoing in the upper cladding 
(Fig. [2)3)7 ^3 = and the proper normalization condi- 
tion is Y3 = 1. All other coefficients can be found from 
transfer-matrix theory, the relevant expressions being ob- 
tained from (ID4|) . (|D5[) by the replacements W — > Y, 
(except in the exponential func- 



X — > Z and gj 
tions). 



Appendix E: MATRIX ELEMENTS FOR 
DIFFRACTION LOSSES 



We calculate the matrix elements ([221) between guided 
and radiation modes of the effective waveguide, to be 
used in formula ([2"Tj) for the loss calculation. The expres- 
sions are similar to those of Appendix [Bj but care must 
be taken to distinguish between the different quantum 
numbers. Guided modes are labelled by the wavevector 
k + G = g and by the index a, which are combined as 
before into a single index /1. Radiation modes are la- 
belled by the wavevector k + G' = g', the frequency u>, 
the polarization A, and an additional index which spec- 
ifies whether the mode is outgoing in medium 1 or 3: 
all these quantum numbers will be grouped into a single 
index r. Like for the matrix elements between guided 
modes, there are four possible combinations of polariza- 
tions. The z integral can be broken into three terms over 
the regions 1,2,3 and it can be expressed in terms of the 
following integrals: 



h± 
h± 



/ e -(*3 (l ±^)(z-d/2) dz = (x3m ± iq3r yi y 
J d/2 



d/2 
d/2 



i(<^±<j2r)z dz = sin((g^ ± q 2r )d/2) 



(g M ± g2r)/2 

J — 00 



d/2 
d/2 



The matrix elements between guided and radiative 
modes are found as follows: 



w ™_te = ^j ^e s -e s ,{(e 1 ) 2 m (G,G')B;, 

X {Wlrh+ + Xirh-) + (e3)^3(G, G')A^ 

x (Wsrh- + X 3r I 3+ ) + {e 2 fr) 2 {G, G') 
x [(A^W 2r + B^X 2r )h- + (A^X 2r + B^W 2r )h+] } 

(El) 

W™"™ = Vi ( G » G')-Di M [{gg + iw^g ■ g')Y lr I 1+ 
+{gg' - ixi»qirg ■ g)Z ir h-\ 

+V3(G, G')C 3 * M [(.g.g' - ix^<iArg ■ g)Y 3r h~ 
+ {99 + *X3m93t-5 • g)Z 3r l3+] 
+r,2(G,G'){(C^Y 2r + D* 2ll Z 2r ) 

x{gg' + Wirg ■ g)h- 

+(C^Z 2r + D* 2li Y 2T ){gg' - q^-g ■ g)h+] , (E2) 
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,i,TE-TM _ ■ / i<> 

ri a r t 



Eg ■ g' {£iVi(G,G')qi r Bt^ 
x(-Y lr h+ + Zxrh-) 

+ e 3 7?3(G, G')q 3r A* 3fl (-Y3rh- + Z 3 rh+) 
+e 2 rj 2 (G,G')q 2r [(-A^Yir + B^Z 2r )h- 
+(A* 2fl Z 2v - B^Y 2r )I 2+ ]} , (E3) 

T.TM-TE w r . „ r _ /„ n » 

?W = — 9 ■ e g ' |-ei»?i(G, G JXim-Dim 
K(Wir/l+ + XirJl-) 

+e 3 ?73(G, G')XSMCa f .(W3r/3- + X 3r I 3+ ) 
-ie 2m {G,G') qii [(D^X 2r - C^W 2r )I 2 - 
+(D* 2ll W 2r - C^X 2r )I 2+ ] } . (E4) 



Appendix F: FOURIER TRANSFORM OF 
TRIANGULAR HOLES 

We give here the complex Fourier transform of the di- 
electric constant in the unit cell for the triangular lattice 
of trian gula r ho les, to be used for the calculations shown 
in Sees. IIIII and IIVI Let us refer to the geometry shown 
in Fig. [5l the side of the triangle and the lattice constant 
being denoted by L and a, respectively. The triangle is 
made of a medium of dielectric constant e\ (in the case 
of an air hole, e% = 1) embedded in a background with 
dielectric constant £2- The origin of coordinates in the 
unit cell is taken to be at the center of the triangle. The 



filling factor of the triangle is 

V3L 2 



f = 



4A ' 



(Fl) 



where A = \/3a 2 /2 is the unit-cell area. The Fourier 
transform is written as 



e(G) = 



e(p)e 



iG p 



dp. 



(F2) 



For zero reciprocal lattice vector it is obviously given 
by: 

e(G = 0)=/ei + (l-/)e a . (F3) 
For nonzero reciprocal lattice vector it is calculated as 

e(G ?6 0) = /(ci - e 2 )[I(G x ,G y ) + I(-G x , G y )], (F4) 

where 

/(G B ,G„ = 0) = 4(l-e- ito )--, (F5) 

gi g? 



I(G x ,G y ^Q) 



— e v 
9y 



-) 



L -iSjL. ( 9x-g v \ . (9x\\ 
(F6) 

with g x = G x L/2, g y = \/3G y L/2. In this work, the tri- 
angles are oriented in such a way that one of the triangle' 
sides is along the x-axis as in Fig. [5] (i.e., along the T-K 
direction of the triangular lattice): if other orientations 
have to be considered, the simplest way is to apply a 2D 
rotation matrix to the vector G before calculating the 
Fourier transform with the above formulas. 
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